Sexual replication in the quasispecies model 
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This paper develops a simplified model for sexual replication within the quasispecies formalism. 
We assume that the genomes of the replicating organisms are two-chromosomed and diploid, and that 
the fitness is determined by the number of chromosomes that are identical to a given master sequence. 
We also assume that there is a cost to sexual replication, given by a characteristic time r see k during 
which haploid cells seek out a mate with which to recombine. If the mating strategy is such that 
only viable haploids can mate, then when T see k = 0, it is possible to show that sexual replication 
will always outcompete asexual replication. However, as T see k increases, sexual replication only 
becomes advantageous at progressively higher mutation rates. Once the time cost for sex reaches 
a critical threshold, the selective advantage for sexual replication disappears entirely. The results 
of this paper suggest that sexual replication is not advantageous in small populations per se, but 
rather in populations with low replication rates. In this regime, the cost for sex is sufficiently low 
that the selective advantage obtained through recombination leads to the dominance of the strategy. 
In fact, at a given replication rate and for a fixed environment volume, sexual replication is selected 
for in high populations because of the reduced time spent finding a reproductive partner. 
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I. INTRODUCTION 

Sexual reproduction is the observed mode of reproduc- 
tion for nearly all multicellular organisms. As such, the 
evolution of sex has been one of the central outstanding 
questions in evolutionary biology. 

One of the biological explanations for the existence of 
sex is that it provides a natural mechanism for diploid or- 
ganisms to eliminate deleterious mutations from a popu- 
lation 1]. The idea is that, by reproducing via a haploid 
intermediate, it is possible for haploids without defective 
genes to recombine with one another, thereby prevent- 
ing the accumulation of deleterious mutations. Other 
explanations that have been advanced are that sex leads 
to greater variability in a population, making the popu- 
lation more adaptable in adverse conditions. It has also 
been postulated that sex evolved as a mechanism for cop- 
ing with parasites p] . 

In recent years, there have been a number of numeri- 
cal studies focusingon the evolutionary dynamics of sex- 
ual replication 0,11 H, 0, ||. These studies have es- 
tablished that, depending on the choice of parameters, 
either sexual or asexual modes of reproduction are the 
advantageous replication strategy. One study in partic- 
ular argues that sexual reproduction is favored when the 
number of daughter genomes produced by the parents is 
high, since this reduces the amount of time required to 
find a reproductive partner 0] ■ 

In this paper, we present a relatively simple evolution- 
ary dynamics model that allows us to compare sexual 
and asexual replication strategies. The model is ana- 
lytically solvable, and treatable within the quasispecies 
formalism . The essential result is that sexual replica- 
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tion is favored in populations with low replication rates, 
and when the characteristic time associated with finding 
a reproductive partner is small compared with the time 
scale associated with replication. These results suggest 
that increasing population density favors the sexual repli- 
cation strategy, since it reduces the time scale associated 
with finding a mate. 

This paper is organized as follows: In the following sec- 
tion (Section II) , we develop a simplified model for sexual 
replication, whose steady-state behavior we proceed to 
characterize in Section III. In Section IV we compare sex- 
ual and asexual replication, and establish regimes where 
each is the preferred mode of reproduction. We conclude 
the paper in Section V with a brief discussion and an 
outline of avenues for future research. 



II. A SIMPLIFIED MODEL FOR SEXUAL 
REPLICATION 

In a simplified model for sexual replication, we as- 
sume that we have a population of single-celled organ- 
isms, where each organism has a genome consisting of two 
chromosomes. We assume that each chromosome may be 
denoted by a linear symbol sequence a = s\ . . . Sl, where 
each letter, or base, Sj, is chosen from an alphabet of size 
S (= 4 for terrestrial life). We further assume that there 
exists a "master" sequence cto for which a given chromo- 
some is functional. It is assumed that a chromosome is 
nonfunctional whenever a ^ <7q (that is, the genes on 
such a chromosome are defective). 

Within this approximation, there are three distinct 
types of genomes in the population: (1) {00, 00} - 
Genomes where both chromosomes are identical to the 
"master" sequence. (2) {(j ,<j ^ er } - Genomes where 
only one of the chromosomes is identical to the "master" 
sequence, while the other chromosome is defective. (3) 



{a ^ ffo.u' ^ do} ~ Genomes where both of the chromo- 
somes are defective. 

We are therefore dealing with a diploid population. If 
we assign the gene sequence rxo & s Viable, while all other 
gene sequences are Unviable, then our three genome types 
may be classified as {V,V}, {V, U}, and {U,U}, where 
V/U stand for Viable/Unviable. 

We assume that the organisms replicate with a first- 
order growth rate constant. For the three distinct 
genome types, the first-order growth rate constants are 
taken to be nyy, kvu, and kjju- We have that nyy > 
kvu > k-uu- 

The sexual replication of the population occurs as fol- 
lows: The diploid organisms divide to form a population 
of haploid organisms. It is assumed that those haploid 
organisms containing a genome of type U are incapable 
of participating further in the reproductive process, so 
that only viable haploids can recombine with each other. 
The newly formed diploids then divide via the normal 
mitotic pathways to form two new daughter cells. 

To develop a set of ordinary differential equations 
governing the replication dynamics described above, we 
first let nyy denote the total number of organisms with 
genome of type {V,!^}, nyu denote the total number 
of organisms with genome of type {V, U}, and nuu de- 
note the total number of organisms with genome of type 
{U, U}. We also let ny denote the population of vi- 
able haploids. We then wish to obtain expressions for 
dnvv/dt, dnvu/dt, dnuu/dt, and dny/dt. 

First note that, the diploid to haploid division leads 
to destruction of each of the diploid genomes at a rate 
given by — Kyynyy for {V, V}, and similarly for the other 
genomes, and a creation of viable haploid genomes at a 
rate given by 2nyynyy + Kyunyu- 

If we let r see fc denote the average amount of time 
a viable haploid spends searching for a viable haploid 
mate, then in a given amount of time dt the total num- 
ber of viable haploids who have recombined is given 
by nydt/r seek (the individual times are Poisson dis- 
tributed). Therefore, recombination leads to a destruc- 
tion rate of haploids given by ny/r seek , and a creation 
rate of diploids given by (l/2)ny /r seek . 

If we let p denote the probability of correctly repli- 
cating a chromosome, then neglecting backmutations we 
have that V — > V with probability p, V — > U with proba- 
bility 1— p, and U — ► U with probability 1. Using this, we 
can construct the various possible replication pathways 
and their associated probabilities, illustrated in Figure 2. 
From these pathways, we can construct the contribution 
to nyy, nyu, and nuu m turn. 

For nyy, all three replication pathways in Figure 1 give 
a contribution. Taking into account probabilities and 
degeneracies, we have a total contribution of 2p 2 from 
the first pathway, 2p(l — p) from the second pathway, 
and (1 — p) 2 /2 from the second pathway, giving a total 
rate of production of {l/2)ny /r seek x [2p 2 + 2p(l — p) + 
(l-pf/2] = {l/4)n v /T seek (l+p) 2 . 

For nyu, similar reasoning gives a rate of produc- 



© (0) 
©— © 

©- 




© 

© 

© 











2p(1-p) 



(1-p)2/2 



(1-p)*/2 



FIG. 1: The various replication pathways and their associ- 
ated probabilities. The factor of 2 in the second pathway 
comes from the fact that either the top or bottom parent "V" 
chromosome can form a daughter "U" chromosome. 



tion of {l/2)ny/r seek x [2p(l - p) + 2(1 - p) 2 /2] = 
(l/2)n V T seek (l -p 2 ). 

For nuu, we obtain (l/2)n v /r seek x (1 - p) 2 /2 = 
(l/4:)ny/T seek (l-p) 2 . 

Putting everything together, we obtain the system of 
differential equations, 
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(1) 



dt T seek 
We can non-dimensionalize these equations by defining 

T = t/T seek , and iiyy — KyyT seek , k VU — KyuTseek, 

iiuu — Kuu T seek- Then we obtain, 

dnyy 



dr 
dnyu 

dr 
dn v 

~dT ~ 



= -Kyynyy + —(1 + p) 

ny 2 
= -Kyunyu + —(1 ~ P ) 

2kyynyy + kyunyu — ny 



(2) 



Define n = nyy +n vu + nuu + ny, n' = n vv +n vu + 
nuu, xyy = nyy/n, x' vv = nyyjn', xyu = n vu /n, 
x'vu = nvu/n', x V u = ™uu/n, x' uv = mu/n', x v = 
n v /n. Note that x' vv — xyy/(l — xy), x' vu — xvu/{^ — 
x v), x'uu = %uu/{l — xy), and n v /n' = x v /(l — xy). 

Then defining k(t) = (l/n)dn/dr, k(r)' = 
(l/n')dn' /dr, we have, 



k(t) = kyyXyy - KuuXuU 

k(t)' = -{kyyx'yy + kyux'yu + kuux'uu) 

Xy 
1 — Xy 



(3) 
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Re-expressing the dynamical equations in terms of the 
xyy, %vu> x uUj x v population fractions, we have, 



dx 



vv 



dr 
dx V u 

dr 
dx vu 



-{kvv + k(t))xvv + -^-(1 +P) 2 
-(kvu + k(t))x V u + -^-(1 - p 2 ) 

Xy 



where k' vv = Kyy/(1 + kyy), Kyu = Rvu/Q + Kvu)- 

We can further simplify the notation by defining k — 
kyy, and a — kyu/kyy. Then R(t — oo)/k is the 
solution to the quadratic, 



= x 2 — A(p, k, a)x — B{p 1 k, a) 



(9) 



- -(%£/ + k(t))xuu + - p) 2 

GET 4 



"S 7 



= -(1 + k(t))xv + 2k vv x vv + k vu x vu 



where, 
A(p, k, a 



(4) 



1 + n v 2 ' 

la 

B p, k, a) = — — — p 

1 + k 1 + an 



' [2(i±^) 2 -l] ' " ^ 



2 1+ ttK 



(1+/) 



(10) 



III. EQUILIBRIUM MEAN FITNESS RESULTS 

At steady-state, the above time derivatives may all be 
set to 0, giving, 

2kyyxyy + kvuxvu 



1 + K(T = OO) 

Therefore, at steady-state, we have that 
1 



(5) 



k(t = oo) 



[((2£yyXyy + kvuxvu) 
Xy 



\—Xy 

-{kyyXyy - kmjXuu)] + ~ 
1 



Xy 



1 — Xy 

R(t = oo) 



Xy + k(t = OO) — Xy(l + k(t — OO))] 

(6) 



Differentiating both sides of the quadratic, it is possible 
to show, after some manipulation, that dx/dp > 0, and 
hence that the mean fitness is an increasing function of 
P- 



IV. COMPARISON OF SEXUAL AND 
ASEXUAL REPLICATION 

If, for simplicity, we assume that Kuu = 0, then 
for asexual replication the steady-state value for R(t — 
oo)/k may be readily characterized 0: It is given by 
max{2( 1 -i £ ) 2 — l,ap}. Therefore, if p CT it is defined by 
the equality 2((1 + p cr it)/2) 2 - 1 = aPcrit, then, 



^(r = oo) _ f 2(1±£) 2 -1 ifpebcrit.l] 
ap if p S (p,Pcrit\ 



(11) 



and so it is equivalent to measure the mean fitness of the 
population at steady-state using either R(t) and R(t)'. 
This is important, because, in comparing mean fitness 
results with an asexually replicating population, the nat- 
ural mean fitness result to use is R(t)'. The equivalence 
between R(t) and R(t)' means that we can compute R(t) 
at steady-state and compare the results directly with the 
value of R(t) for the asexually replicating population. 

Plugging the steady-state value of xy into the steady- 
state equations for xyy and xyjj, we obtain, 



We can compare these values with that obtained for 
sexual replication. 



A. Case 1: k = 

We begin by considering the case where there is no 
time cost associated with sex, so that r see k — =>■ R = 0. 
Then 



•1+^2 



1 



o = [r vv (2(— ±y 

2 1 + K(T = OO) 

, l/i , \2 K VU X VU 

+ 4 {1+P) l + ft(T = oo) 

= [*M~(1 -P 2 ) T— ~T \ 

2 1 + K(T = OO) 



+ (1-P 2 )- 



1 



-KyyXyy 



1) — k(t = (X))]xyy 

1) — R(t = oo)]xyu 
(7) 



A(p,0,a) = [2(I±Z) 2 -l] 
B(p, 0, a) = ap 



1 



a(l+p 2 ) 



(12) 



1 + R(t) 

We can solve the first equation for xyu in terms of xyy 
Plugging the resulting expression into the second equa- 
tion, we obtain, after some algebra, the quadratic, 

= R{t = oo) 2 

2 2 

(8) 



I^y(2(^) 2 - 1) - ^k' vu (l+p 2 )}R(r = oo) 



kyykyjjP 



We claim that for k = 0, k(t = oo)/« > 
max{2( -!-ip) 2 — 1, ap}, with equality occurring only when 
p = 1 for arbitrary a, and a = 0, 1. 

To prove this claim, note first that for p = 1, we have 
2(-Mp) 2 — 1 = 1, and that R(t = oo)/n = 1 = max{l, a}. 
So since the claim is true for p — 1 , we may now consider 

pe[o,i). 

If a = 0, then k(t = oo) = max{2(±±£) 2 - 1, 0}, while 
if a = 1, then «(r = oo) = p = max{2(-!-rp) 2 — l,p}, since 
2(ii£) 2 — 1 < p for p e [0, 1], with equality occuring only 
when p = 1. 

So, we now consider the case where a € (0,1), and 
p G [0,1). If we define k(p) = 2(±±£) 2 - 1, then we 
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have two possibilities: Either max{fc(p), ap} — k(p), or 
max{fc(p), ap} = ap. We will consider each of these two 
cases in turn. 

So, first assume that max{fc(p), ap} — k(p). We wish 
to show that, 



\[k{p)~a{l+p 2 )+yj (k{p) - +p 2 )f +4ap] > k(p) 

(13) 

After some manipulation, we obtain that this condition 
is equivalent to the condition that 



> / + 2p 3 - 2p - 1 



(14) 



To establish this inequality for p £ [0, 1), note that p 4 + 
2p 3 - 2p - 1 = (p - l)(p 3 + 3p 2 +3p+ 1). Since p 3 + 
3p 2 + 3p + 1 > for p e [0, 1), and since p — 1 < for 
p £ [0, 1), the inequality follows. 

So now suppose that ma,x{k(p) , ap} = ap, so that 
k(p) < ap. Then for our calculations, we first re-write 
A(p, 0, a) as a(p — 1) + (1 — a)k(p). Then we wish to 
show that, 



1 



*(p- 1) + (1 - a)k(p) + 



> ap 



1) + (1 - a)k(p)) 2 + lap] 



(15) 



After some manipulation, this becomes equivalent to the 
condition that, 



k(p) + 1 > a(k(p) + 1) 



(16) 



which is certainly true, since k(p) + 1 > for p e [0, 1), 
and a € (0, 1) by assumption. 

Therefore, we have proven our claim, and hence, within 
this model, sexual replication leads to a greater mean fit- 
ness for a population than asexual replication, assuming 
that there is no cost associated with sex. Figure 2 shows 
a plot of R(t = oo)/k for both sexual and asexual repli- 
cation, assuming a = 1/2 and k = 0. 

Note that for k = 0, if a sexually and asexually repli- 
cating population were placed in an identical flask, then 
under the circumstances dictated by our model the sex- 
ually replicating population would eventually dominate 
the population (that is, the fraction of sexually replicat- 
ing organisms would increase to 1, while the fraction of 
asexually replicating organisms would decrease to 0). 



Case 2: n > 



have, for sexual replication, that 
for p = 1. Since for asexual repli- 



For k > 0, we 
k(t = oo)/k = j±- k 
cation we get k(t = oo)/k = 1 for p = 1, it follows by 
continuity that there exists a regime [p=(n), 1] for which 
asexual reproduction leads to a greater mean fitness of 
the population than sexual reproduction. Presumably, 
as k increases, p={n) should decrease. 




FIG. 2: Comparison of k(t — <x)/k for both sexual and asex- 
ual replication, with k = and a — 1/2. Note that sexual 
replication outcompetes asexual replication for all mutation 
regimes. 



We can determine, for a given k, the mutation regime 
where asexual replication outcompetes sexual replication, 
and the mutation regime where sexual replication out- 
competes asexual replication. To do this, we do not at- 
tempt to compute p=(k) directly. Rather, as a function of 
p, we seek to determine n={p), the value of k for which 
asexual and sexual replication yield identical mean fit- 
nesses. Since by definition k = (p(k)) = k, the function 
K = (p) may be inverted to obtain p = {n). In what follows, 
we restrict our analysis to the case where a £ (0, 1) (if 
a = 0,1, then for k > asexual reproduction outcom- 
petes sexual reproduction for all values of p) . 

Because the mean fitness of an asexually replicating 
population falls into two distinct regimes defined by the 
cutoff Pcrit (at least, within the context of our model), 
we must determine n = (p) separately for p < p crit and 
P > Pcrit- For p < p cr it 7 we have for asexual replication 
that k(t = oo)/k = ap, and hence we must have, 







i 1 

-(ap) 



+ an = 



(l+p 2 ))(ap) 



1 + k = 1 + an. 



(17) 



Assuming that ap > 0, this expression may be re- 
arranged and simplified to, 







{a K= ) 2 p + {ap + 1)(ok=) --(!-<*)(!+ pf (18) 



so that an = , and hence k = , may be solved using the 
quadratic formula. 

We claim that k = is an increasing function of p on 
the interval [0,p cr it]. We can prove this by showing that 
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an = is an increasing function of p on the interval [0,p cr ii]. 
Defining x(p) = an = (p), we have, 

= px(pf + (ap + l)x(p) - ~(1 - a)(l +pf (19) 

Differentiating with respect to p, we obtain, 

= (2px{p)+ap+l)x'{p) + x(p) 2 + ax(p)-(l-a)(l+p) 

(20) 

so we wish to show that x(p) 2 +ax(p) — (1 — a)(l+p) < 
for p £ (0,p cr ;t). Multiplying both sides of the inequality 
by p, and noting that px{p) 2 = \(1 — a){l + p) 2 ~ {ap + 
l)x(p), we have that for p > we need to establish the 
inequality, 




PmCIM) Pm(1'2> 



> -(l- a )(l-p 2 ) 



(21) 



FIG. 3: k= versus p for a = 1/4, 1/2,3/4. A graph of K cri t(a) 
is included as well. 



To prove this inequality, note first that x(0) = a), 
and x'(0) = ~(1 — a)(3 — a) > 0. By continuity, x'(p) > 
in a neighborhood of p = 0. If a;'(p) < for some 
p > 0, then by the Intermediate Value Theorem there 
exists at least one p > for which a;'(p) = 0. If p* = 
inf{p £ [0,l]\x'(p) — 0}, then by continuity it follows 
that x'(p*) = 0, and hence p* > 0. Therefore, x'(p) > 
for p e [0,p*), otherwise by the Intermediate Value 
Theorem there would exist ap** < p* such that x'(p**) — 
0, contradicting the definition of p* . But, since x'(p) > 
for p £ [0,p*), it follows that x(p) is increasing on [0,p*], 
hence x(p*) > x(0) = |(1 - a) > |(1 - a)(l - (p*) 2 ), 
which implies that x'(p*) > =><=. Therefore, x'(p) > 
for p S [0,1], hence on [0,p cr it], n=(p) increases from 
|(1 - a)/a to «;=(p C rit(a))- 

Now, for p € \pcriti 1]) w e have for an asexually repli- 
cating population that k(t = 00)/ n — k(p), hence, in 
this regime, k=(j>) is defined by, 



= k(p) 2 - (- 



1 



- K = 

a 



-k(p) - 



1 



2 1 + an = 



1 + K— 1 + aK = 



(l+p 2 ))fc(p) 



(22) 



which after some manipulation may be re-arranged to 



= k 2 - -\- B(p, a)n = — C(p, a) 



(23) 



where B(p,a) = 1 + l/a + (1 — p)/k{p), C(p,a) 
i^y(l + l/fc(p))(l -p 2 ). We then have that, 



/s=(p)=5(Wl + 4 



'S 2 



1) 



(24) 



We claim that k = (p) is a decreasing function of p for 
P & [Peru, 1]- We will prove this by showing that B and 
C/B 2 are both decreasing functions of p for p £ [p cr it, 1]- 

To prove that _B is a decreasing function of p for p £ 
[Pcrit, 1], we need show that (1 — p)/k(p) is a decreasing 



function of p for p £ \p cr u , 1] • Differentiating, we obtain, 



d l-i 
dp fc(p) 



1 1 
2 _ 



- p 2 + 2p 
k(p) 2 



< 



(25) 



so _B is certainly a decreasing function of p for p S 
[Pcrtt, !]• 

Now, after some manipulation, we can show that, 



C_ 

S 2 



1 2 

—a 

2 



l-p 2 



( l + fc( p) )( l + a __i+^)2 



(26) 



Note that 1 — p 2 is decreasing for p £ [p crit ,l], and 
that 1 + fc(p) is increasing. We also have that, 



d , 1 + ap ^ 
cip 1 + k(p) 



(a/2)p 2 + p+ (1 - a/2) 
(l + fc(p)) 2 



< (27) 



so that (1 + ap)/(l + fc(p)) is a decreasing function of p. 
Therefore, 1 + a — (1 + ap)/(l + k(p)) is an increasing 
function of p, hence C/B 2 is decreasing for p S [p C Ht, 1], 
as we wished to show. 

Figure 3 illustrates the behavior of k- (p) for three val- 
ues of a. 

In what follows, we shall change our notation slightly 
to explicitly indicate that n = also depends on a. Thus, 
we shall re-denote «=(p) by K=(p, a). This notation was 
not needed in the previous arguments, since we were con- 
sidering the behavior of n = at a fixed a. 

We may now summarize the behavior of k = (p, a) as a 
function of p: 

From to p cr it, ^=(p, ct) increases from ^(1 — a)/a to 
Kcrit(a) = K=(Pcrit(ct),at), while from p crit to 1, K = (p, a) 
decreases from K cr it(ct) to 0. This behavior leads to three 
distinct regimes of k. 

For k £ [0, i(l — a) /at], there exists only one value of 
p for which asexual and sexual replication yield identical 
mean fitness results. This value of p is contained in the 
interval [p cr it: !]• As k increases from to ^(1 — a) /a, 
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this value of p decreases. For these values of k, asexual 
replication is advantageous over sexual replication at low 
mutation rates. However, there is a crossover replica- 
tion fidelity where sexual replication becomes advanta- 
geous. As k increases, this crossover replication fidelity 
gets pushed to lower values. This makes sense, since a 
higher value of k corresponds to a greater penalty asso- 
ciated with sex. 

For k G [i(l — a) /a, K cr i t (a)], there exist exactly two 
values of p for which asexual and sexual replication yield 
identical mean fitness results. One value of p is contained 
in the interval [0,p cr it], while the other value is contained 
in the interval \p cr it, 1]- As k increases from ^(l — a) /a 
to K cr it(a), the value of p in [0,p cr it] increases from 
to p crit (a), while the value of p in [p cr jt,l] decreases 
to Pcrit(ct). For these values of k, asexual replication 
is also advantageous over sexual replication at low mu- 
tation rates. As with the previous regime, there is a 
crossover replication fidelity where sexual replication be- 
comes advantageous. However, in contrast to the first 
k regime, there is a second crossover replication fidelity 
where asexual replication again becomes advantageous. 
For these values of k, the cost associated with sex is still 
sufficiently low that sexual replication can become the ad- 
vantageous strategy at higher mutation rates. However, 
the cost of sex is sufficiently high that, at even higher mu- 
tation rates, sexual recombination no longer offsets the 
production of unviable chromosomes from viable ones to 
an extent that makes the strategy advantageous. 

Finally, for k e (K crit (a) , oo) , the cost associated with 
sex is so high that sexual replication is never the advan- 
tageous strategy. 

As a final note for this subsection, we can show that 
Kcrit {(*) is a decreasing function of a. Differentiating the 
quadratic equation given by Eq. (19) with respect to a 
(where x = x(a) = an crit (a)), we have, 

dx 

= (2x Pcrit + ap cr i t + 1) — 

da 

(x 2 + ax - (1 - a)(l +Pcrit))— 

da 

+PcritX + - (1 + Peril) 1 (28) 

Differentiating both sides of the equality defining p cr u , it 
is possible to show that dp cr i t /da — p cr it/{^- — & + Peru)- 
We also have, from a previous analysis, that px 2 + axp — 
(1 — a)p{l +p) = \{1 — a)(l — p 2 ) — x, so that, to show 
dx/da < 0, we need to show that 

i(l - a)(l — p 2 it ) — x 1 x9 

< 2V - c " tJ +p cri tx+-{l+ Pc „t) 2 (29) 

1 - a + p cr i t 2 

Multiplying by 1 — a+p crit and simplifying, this is equiv- 
alent to the inequality, 

< (l-a-x)+p C rit[(l-Ot)(l+x)+X Pc rit + ^(l+Pcrit) 2 } 

(30) 



Since we showed that the expression for x(p) in Eq. (19), 
valid over p G [0,p C rit], is an increasing function of p 
for p e [0, 1], then solving the quadratic in Eq. (19) for 
p = 1 gives x < 1 — a. Therefore, 1 — a — x > 0, hence 
the inequality holds. 

We have therefore shown that an crit (a) is a decreasing 
function of a, hence K crit (a) is a decreasing function of 
a. When a = 1, p cr it = 1, so k = = 0. As a — > 0, 
p cn t -> y/2- 1, so at -» (VV2-3- 1)/[2(>/2 - 1)] 

K = — ► OO. 



C. Consideration of nuu > 

When Kuu > 0, the results for sexual replication re- 
main unchanged. However, the results for asexual repli- 
cation change somewhat, since now an additional lo- 
calization to derealization transition can occur, once 
Kuu = max{Ky V (2((l +p) 2 /2) - l),K VU p,Kuu}- It is 
therefore possible to have a situation where sexual repli- 
cation only becomes advantageous once the mean fitness 
of an ascxually replicating population is Kuu- With an 
appropriate choice of parameters, it is then possible that 
the mean fitness of the sexually replicating population is 
less than Kuu , so that sexual replication will never be the 
preferred replication strategy. We leave the investigation 
of this phenomenon for future work. 

In any event, for Kuu > 0, asexual replication will 
become the advantageous mode of replication at suffi- 
ciently high mutation rates, since the mean fitness of the 
sexually replicating population decreases to zero, while 
after complete derealization over the genome space has 
occurred, the mean fitness of the asexually replicating 
population becomes Kuu- This result, however, is likely 
due to a mating strategy that essentially "throws away" 
the "U" chromosomes. Other mating strategies, where 
all haploids are capable of mating, may not exhibit the 
same effect. 



V. DISCUSSION, CONCLUSIONS, AND 
FUTURE RESEARCH DIRECTIONS 

This paper developed a simplified model for sexual 
replication, and showed, within the context of the model, 
that a sexually replicating population will outcompete an 
asexually replicating one when there there is no cost as- 
sociated with sex. We further showed that if the cost 
associated with sex is sufficiently low, then sexual repli- 
cation becomes advantageous at higher mutation rates, 
because recombination prevents the accumulation of de- 
fective mutations in the diploid genomes (assuming that 
there is a fitness penalty associated with the defective 
mutations). The cost for sex was measured by the di- 
mensionless parameter k, defined to be the product of 
the first-order growth rate constant of the mutation- 
free genomes (kvv), and the characteristic time associ- 
ated with finding a recombination partner (r see fc). Since 
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Kyv — 1/T~rep, where r rep denotes a characteristic repli- 
cation time, it follows that k = T see ^/T rep . Therefore, 
the cost associated with sex is measured by the ratio of 
the time a haploid spends finding a recombination part- 
ner with the time scale for replication. The smaller this 
ratio, the smaller the fitness penalty incurred by repro- 
ducing via a haploid intermediate, and the greater the 
selective advantage for sex. 

The implications of this model are that sexual replica- 
tion is favored in environments where organisms replicate 
relatively slowly, and in environments where the time 
spent finding a recombination partner is small compared 
with the time scale for replication. Thus, sexual repli- 
cation is favored in environments with high population 
density. These results are therefore consistent with the 
observation that sexual replication is the preferred (and 
generally the only) mode of reproduction for nearly all 
multicellular organisms. 

That sexual replication only becomes the preferred 
mode of reproduction at low k suggests why sexual repli- 
cation occurs as a stress response in some organisms, such 
as Saccharomyces cerevisiae (Baker's yeast). When con- 
ditions are favorable, KvVi & n d hence n, are relatively 
high, so asexual replication is the advantageous strategy. 
Under sufficiently adverse conditions, n can drop to lev- 
els where the sexual strategy becomes advantageous. The 
rcplicative strategy that can adopt the optimal replica- 
tion strategy for the given environment will have a selec- 
tive advantage (assuming that resource costs for main- 
taining this switching behavior are not prohibitive), and 
so organisms carrying this strategy in their genomes will 
dominate the population. 

However, as one moves toward more complex life forms, 
the replication rate drops to values such that asexual 
replication is almost never the preferred mode of repro- 
duction, so that the ability for an organism to switch be- 
tween the two modes of reproduction disappears. At this 
point, we postulate that the division of haploid cells into 
two distinct types of gametes, and then later the division 
of the organisms themselves into male and female, are 
the result of selection for evolutionary pathways leading 
to the division of labor and specialization of tasks associ- 
ated with sexual replication. When replication rates are 



low, and when the time cost associated with sex is low, 
then it is likely more efficient (in terms of resource uti- 
lization) to divide the reproductive tasks associated with 
sexual replication among two types of organisms ( "male" 
and "female" ) . The relative fitness advantage as a result 
of such savings in resource costs likely increases with the 
complexity of the organism, leading to a stronger selec- 
tion pressure for a male-female split as organismal com- 
plexity grows. 

In this paper, we assumed that only haploid cells with 
viable chromosomes are capable of engaging in sexual re- 
combination. This allowed a simplified analysis within 
the standard quasispecies formalism. While we obtained 
a selective advantage for sexual replication using this 
mating strategy, a fuller analysis will require the con- 
sideration of various mating strategies on the selective 
advantage for sex. An important such mating strategy, 
which is the opposite of the one considered in this paper, 
is the random mating strategy, whereby all haploids are 
capable of engaging in sexual recombination, and do so 
with a pairwise distribution given by the Hardy- Weinberg 
equilibrium. 

In this vein, one interesting question is to determine, 
for a given fitness landscape, whether there always ex- 
ists a mating strategy for which sexual replication will 
outcompete asexual replication. Additionally, while this 
paper implicitly assumed that the strategy for sexual or 
asexual replication is inherited, future studies should con- 
sider genomes where genes for sex are explicitly included. 
This leads to the ability for sexual organisms to mutate 
into asexual ones. As the selective advantage for sexual 
replication disappears (as a function of mutation rate, 
for instance), the models may exhibit localization to de- 
localization transitions over the portions of the genome 
controlling sex. 
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